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ABSTRACT 

This invention is a method of separating induced microseismicity signals from seismic 
signals acquired within active seismic monitoring operations carried out in underground 
zones under development which has application for underground hydrocarbon or fluid 
storage reservoir. - Seismic records are formed from signals emitted by one or more 
seismic sources controlled by orthogonal signals. In this case, the signals are processed to 
separate the respective contributions of the seismic sources to the signals received and to 
reconstruct the seismograms equivalent to those that would be obtained by actuating the 
seismic sources separately. The induced microseismicity signals (passive seismic) are 
separated from the seismic signals coming from sources (active seismic) essentially by 
isolating the contribution thereof by comparison with a spectral reference model at the 
ftmdamental frequencies emitted and at the respective harmonics thereof, and by 
reconstructing in the time domain the induced microseismicity signals and the signals 
coming from seismic sources. - Applications : seismic monitoring of underground 
hydrocarbon or fluid storage reservoirs. 
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FIELD OF THE INVENTION 



The invention relates to a method intended for discrimination of microseismic 
signals induced by the development of an undergroxmd reservoir, from among seismic 
signals emitted within the context of seismic monitoring operations. 



is well-known to permanently install seismic sources and seismic receivers, at the 
sur&ce or in wells, whether reservoir e3q)loration or development wells, and to carry 
out, at predetermined time intervals, recording of the seismic signals from the 
formation. These signals may be seismic signals induced by the site development 
10 activity or reflected by the geologic discontinuities of the formation in response to the 
emission of disturbances in the ground by means of one or more seismic sources. 
Comparison of the recordings performed at different times of the site development 
provides information on the evolution thereof. 



15 Various long-term seismic monitoring systems are described for example in patents 

EP-591,037 (US-5,461,594), FR-2,593,292 (US-4,775,009), FR-2,728,973 (US- 
5,724,3 11) or FR-2,775,349. 

Patent EP-748,457 (US-5,724,311) filed by two of the applicants describes a 
permanent monitoring system allowing to obtain perfect reproducibility of the operating 
20 conditions in repetitive seismic monitoring operations carried out in an underground 
zone traversed by at least one well or borehole, notably in an undergroimd gas storage 
reservoir. This system conq>rises, installed on a permanent basis, one or more set(s) of 
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In order to follow the evolution of an undergroimd reservoir under development, it 



BACKGROUND OF THE INVENTION 
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seismic receivers (buried at the sur&ce or placed in one or more wells), several 
repetitive seismic sources (buried or at the surfece), and a permanent linking network 
for selective energy supply to these sources. A central station remote controls 
selectively each one of these sources and records the seismic signals coming from the 
5 underground zone in response to the seismic waves transmitted selectively in the groimd 
by the sources. 

All these permanently installed soiirces whose coupling with the surrounding 
formations remains stable and this supply network, at least partly buried and whose 
sur&ce coverage area is limited, allow to carry out a series of long-term seismic 
10 monitoring operations imder stable operating conditions, without any risk of 
incompatibility with the development site activities. 

Patent FR-2,775,349 (US-6, 182,0 12) filed by the same two applicants describes a 
method and a device intended for permanent monitoring of an underground formation in 
- which cavities are formed for seismic sources (preferably deep enough to be coupled 

15 with the formation below the weathered zone (WZ)), at least one well is drilled either 
below each source or in the immediate vicinity thereof An antenna consisting of a 
plurality of seismic receivers (geophones, hydrophone, etc.) is lowered into each well 
thus drilled and coupled with the surrounding formation. With this device, emission- 
reception cycles are carried out with emission of seismic waves in the formation by at 

20 least one source and acquisition of the signals picked up by the seismic receivers, in 
response to the waves reflected by the formation. The receivers can for example be 
associated with a fluid transfer tube lowered in the well in order to connect an 
underground zone to a surface development device and coupled with the formation 
surrounding the well. - 
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The various sources of the seismic monitoring device can be actuated successively 
with a sufficient time mterval between the triggering times for reception of the waves 
reflected by the investigated zone. It is also possible to use several seismic sources 
emitting the same signals that are triggered simultaneously to increase the power 
5 emitted. 

Patent FR-2,805,051 filed by the applicants describes a method and a system 
. intended for seismic monitoring of an underground formation. As diagrammatically 
shown in Figs. 1 to 4, the system comprises for example a network of seismic antennas 2 
consisting each of a series of seismic pickups 4 arranged at regular intervals along a 

10 well 3 drilled in the ground. This network can be regular as shown in Fig.2 or irregular. 
The pickups are, for example, vertically oriented one-directional or multi-axis 
(triphones) geophones and/or hydrophones. A seismic source S is arranged near to each 
antenna 2. Piezoelectric type vibrators such as those described in patent FR-2,791,780 
(US-6,338,394) filed in the joint names of the applicants are advantageously used as 

15 sources and permanently installed in the immediate vicinity of each antenna 2. An 
electronic control and recording module 6 (Fig.l) can be associated with each anteima. 
These electronic modules 6 can be connected to a central control and synchronization 
laboratory 8. According to another implementation mode diagrammatically shown in 
Fig.4, antennas 2 are directly connected to a central laboratory 8 including all the 

20 individual electronic control and recording modules 6 (Fig.l). 

Seismic waves are emitted in the underground formation by one or more seismic 
sources (vibrators) 5. In the case of multiple sources, the latter emit simultaneously; 
they are controlled by orthogonal signals so as to form a composite vibrational signal. 
The vibrators work permanently or at least for a suflBciently long time to x>bta]n' a 
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sufiQcient signal-to-noise ratio. The seismic waves generated by seismic sources 5 are 
propagated downwards (downgoing waves 9). These incident waves are first recorded 
by receivers 4 of each well 3. The waves reflected by the discontinuities of the zone 
(seismic interfaces) are propagated upwards. These upgoing waves 10 are also recorded 

5 by the various receivers 4. The upgoing and downgoing waves are thus superimposed in 
the seismograms. The records are subjected to the conventional processings known to 
the man skilled in the art suited for those obtained by means of the VSP technique 
(Vertical Seismic Profiling). Finally, the respective contributions of the seismic sources 
to the composite vibrational signal are discriminated and the seismograms equivalent to 

10 those that would be obtained by actuating the seismic sources separately are 
reconstructed. Sinusoidal signals of different fi-equencies, in their fimdamental 
components as well as in their respective harmonics, or signals based on wavelets, 
Legendre polynomials or random series, etc., are for example used as orthogonal 
signals. 

15 In the case notably where the orthogonal signals emitted are sinusoids, 

discrimination of the respective contributions of the seismic sources is for example 
performed by detemoining the amplitude and the phase of the conposite vibrational 
signal at the fundamental fi-equencies of the control signals applied to the seismic 
sources. This discrimination comprises for example weighting of the signals recorded 

20 by a bell weighting (or apodizing) factor and determination of the amplitude and of the 
phase of the composite signal, or selection by Fourier transform of lines of the complex 
spectrum associated respectively with the various weighted signals. Reconstruction of 
the seismograms corresponding specifically to the various seismic sources is carried out 
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for example by applying, after separation thereof an inverse Fourier transform to the 
lines associated respectively with the various weighted signals. 

Furthermore, it is well-known that reservoir development induces mechanical stress 
variations in the reservoir and in its overlying bed. These stress variations have the 
5 efifect of generating seismic signals that are propagated in the formations. The receivers 
of the various antennas 2 record these signals (upgoing waves 11). 

Since seismic sources S are arranged much closer to receivers 4 of antennas 2 than 
to the reservoir and since the energy they emit is much more intense than the energy of 
the induced microseismic signals generated in the reservoir, the microseismic signals 
10 are generally concealed and imdetectable. 

SUMMARY OF THE INVENTION 

The method according to the invention allows to separate induced microseismicity 
signals and multisource seismic signals within the context, of monitoring operations in 
an underground zone imder development. 

15 It comprises emission of seismic waves in the formation by coupling therewith one 

or more seismic sources. In the latter case, these sources emit simultaneously 
orthogonal signals so as to form a composite vibrational signal, reception of the signals 
reflected by the formation in response to the emission of seismic waves, recording of 
the signals received by at least one seismic pickup and processing of the signals 

20 recorded so as to separate the respective contributions of the seismic sources to the 
signals received and to reconstruct the seismograms equivalent to those that would be 
obtained by actuating the seismic sources separately. 
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The induced microseismicity signals (passive seismic) are separated from the 
seismic signals emitted essentially by isolating the contribution of each source by 
conq)arison with a reference spectral model at the fundamental frequencies emitted and 
at the respective harmonics thereof and by reconstructing in the time domain the 
5 induced microseismicity signals and the seismic signals emitted. 

The spectral contribution of the microseismicity signals to the spectrum of the 
signals received is obtained for example by subtracting the anq)litude and phase values 
associated with the reference spectral model from the amplitude and the phase of the 
spectrum associated with the records. 

10 According to an ini^lementation mode, a current model formed by updating a 

spectral model taking account of the contribution of previous recording cycles is for 
example used as the reference spectral model. 

The current spectral model can be formed by determining a mean value of the 
frequency spectra formed from (earlier and/or later) records obtained for the same 
15 source and the same frequencies. 

According to another implementation mode, the current spectral model is formed by 
determining a median value of the frequency spectra formed from earlier records 
obtained for the same source and the same frequencies. 

According to another implementation mode, the current spectral model is formed by 
20 extrapolation or interpolation from the frequency spectrum, from close spectral values. 

According to another implementation mode, the current spectral model is formed by 
using the line obtained at this frequency. 
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According to an implementation mode, discrimination of induced microseismicity 
signals from among signals emitted within the context of active seismic monitoring 
operations is obtained by means of the following stages : 

a) for each recording n of a recording cycle p, the respective contributions of the 
5 various sources at the fundamental frequencies are calculated, 

b) the ratio of the contribution to a current spectral model formed by updating a 
previous spectral model from frequencies emitted during the previous recording and 
from their harmonics is then calculated, 

c) the part of the recording n of cycle p that can be associated with the active 
10 seismic monitoring operations is deduced, 

d) the part of the recording n of cycle p strictly relevant to the passive microseismic 
activity is deduced, 

e) the seismograms that can be associated with the active seismic monitoring 
operations are formed by inversion in the time domain of the respective spectral 

15 contributions of each seismic source at the fimdamental frequencies and at their 
harmonics, and 

f) the underlying microseismic signals contained in the records are formed by 
inversion in the time domain from the part relevant to the passive microseismic activity. 

The spectral contribution is obtained by multiplying a transfer fiinction between a 
20 wavelet characteristic of the source and a seismogram associated with a given receiver 
by a wavelet characteristic of the source. This transfer fimction can be continuously 
updated during a current cycle from an estimation made during a previous cycle and 
from an initial estimation made during the current cycle by the relation. 
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BRffiF DESCRIPTION OF THE FIGURES 

Other features and advantages of the method according to the invention will be 
clear from reading the description hereafter of an embodiment given by way of non 
limitative example, with reference to the accompanying figures wherein : 

5 - Figure 1 diagraninmtically shows the monitoring device installed in the field for 
seismic monitoring of an underground formation, 

- Figure 2 shows a mode of layout of the emission-reception units in the field, 

- Figure 3 diagrammatically shows one of these emission-reception imits comprising a 
seismic source and a reception antenna, 

10 - Figure 4 shows a variant of the monitoring device of Fig. 1 , and 

- Figure S shows examples of signals and frequency spectra mentioned in the 
description. 

DETAILED DESCRIPTION 

The method thus allows to carry out seismic monitoring operations in an 
15 imderground zone by means of a series of seismic pickups and one or more vibrational 
seismic sources. If several simultaneous sources are used, they are actuated 
simuttaheously by signals at different frequencies selected so' as to allow discrimination 
of the contributions of each source in the seismograms obtained from the signals 
received and recorded. 

20 This is generally performed by control of the various sources by sinusoids of 

different frequencies and by means of known numerical calculation techniques such as 
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inverse Fourier transform to separate the contributions to the seismograms obtained 
from the various seismic sources. 

Distortion cannot be disregarded with mechanical sources. Simultaneously with a 
frequency fj, a single source S, emits frequencies 2^, 3^ ... nfj. Consequently, if 5 and ^ 
5 are the respective frequencies of the two sources Si and Sj of the series of sources, we 
must have fi 9^ ^, as well as ^ 2fj, fj 5^ 3£j, . . . fi n^ in order to be able to separate their 
respective contributions. 

As already described in the aforementioned patent FR-2,805,051, if the con:q)osite 
signal Pt consisting of the sxmi of N sinusoids Aj, Oj} with 1 < i < N is emitted, all 

10 the frequencies being contained in a spectral band contained between two boundary 
frequencies and f^ the seismogram Tt observed at point R will have as Fourier 
transform at the frequency fi the number of module Aj and of phase Oi equal to the 
amplitude and to the phase of sinusoid Tj. It is thus possible to reconstruct by iuverse 
Fourier transform seismogram Tt by emitting successively all the sinusoids of 

15 frequencies 5> to ff. 

In the case where, for example, all the amplitudes Aj are 1 and all the phases are 
zero, the signal Pt obtained is very close to the signal resulting from autocorrelation of a 
sliding-frequency signal contained in the sweep range [ft - jQ, commonly used in 
vibroseismic prospecting. According to the discrete Fourier transform theory well- 
20 known to the man skilled in the art, if one wants to Usten to source Si during the time te, 
the frequency increment between the sinusoids is Af = l/t© and the number of sinusoids 
necessary is Nf= (ff- ft)te. 
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N vibrators SI, S2, Si,... Sn installed in the field can thus be simultaneously excited 
by means of vibrational signals whose fi-equencies are such that each source is excited 
successively by each above sinusoid Nf at any time, on condition that the respective 
frequencies of the sinusoids emitted at the same time by the various seismic sources are 
5 all different from one another. Nf recordings corresponding to the Nf frequencies 
selected in the sweep range are thus performed with each vibrator. Separation of the 
signals received by the pickups in the field, in response to the simultaneous emission of 
the various signals, is thus obtained by selection of the line at the suitable frequency. 

Cycles of Nf recordings are thus repeated. The record En of a given acquisition cycle 
10 contains a set of ftindamental and harmonic frequencies which only depends on n. 
StartiQg from this observation, the method that will be applied to detect the effects of 
microseismicity is based on the observation which has proved tme in practice according 
to which the amplitude and phase variations observed at the frequency f and at the tune 
t are very httle different from those observed at the frequency f - Af and at the time t - 
15 At. This allows to form a cumulative spectral model constructed fit)m records of 
previous cycles at the same frequencies but formed on a long time during which the 
characteristics of the seismic signals emitted are likely to vary, and permanently adapted 
from records obtained during a sufficiently short period for the variations of the signal 
to be disregarded. The efficiency of the proposed method is optimum when the 
20 successive frequencies emitted by the same source are regularly incremented. Cycles in 
upgoing and downgoing frequencies can for example be alternated. 

The procedure for showing in practice this spectral contribution of the vmderlying 
microseismic signals is as follows. 
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Notations : 

■ fi^ designates the fundamental frequency emitted by source i during the recording 
n (this frequency is independent of the cycle). 

In order to simplify the notations of the following quantities, the subscript giving 
5 the number of the receiver is left out. Whatever this number, the quantities are 
calculated in the same way. 

■ Rp^ designates the recording of rank n of cycle p, 

■ Cp^^n the contribution of source i to the recording n of cycle p, 

- Mp^ the initial model of the active part of the recording n of cycle p. This model 
10 contains only the frequencies emitted and their harmonics during recording n, 

• Ep^ the ratio between the active contributions of the recording n of cycle p and the 
initial model of the active part of the recording n of cycle p, 

- Ap^ the active part of the recording n of cycle p, 

■ Pp^ the passive part of tiie recording n of cycle p (microseismic), 

15 ■ k the differences updating coefficient (typically of the order of 1 5 to 25 %), 

- h the models updating coefficient (typically of the order of 5 to 10 %). 

Considermg these notations, the method first consists, for each recording n of cycle 
p, in calculating the respective contributions Cp,u of the various sources at their 
fundamental frequencies and at their harmonics by means of the method described in 
20 the aforementioned patent FR-2,805,051. Owing to the possible presence of 



12 



microseismic energy, the contribution of the various sources is likely to be affected by 
noise. 

This estimation can be improved by associating the contributions of all the receivers 
of a same recording Cp,i,„. These contributions can be considered to be the product of 
two &ctors among which one depends only on the source (and not on the receiver) and 
the other on the receiver, and not on the source. 

If subscript r represents the receiver, 

T only depends on p, i and n by the variation of the characteristics of the receiver 
and of its environment with time. These variations are very slow and can be disregarded 
or estimated more precisely. 

S and T are determined to within one &ctor. The average of the downgoing waves 
can thus be selected for Spj^. This average is obtained by means of a calculation well- 
known to those skilled in VSP processing (flattening of the downgoing waves and 
average). The Tp.i^^ are then the transfer fimctions between the downgoing wave and the 
seismogram. 

In cases where the variation of T is to be taken into accoimt, updating coeflBcient h 
can for example be used to update Tp^i^r from Tpj^.i.r as foDows : 

If Cop,!^ designates the initial estimation of contribution C and TOp i n ^ that of the 
receiver factor, S is calculated by flattening and averaging of the downgoing waves, 
then TOp^i^ is calculated by the relation : 




Co 




The final estimation of the receiver fector is obtained by : 



Tp,i,n,r — (1 - h)Tp,i^.i^ 



+ hTo, 



and the final contribution of source i to cycle p of recording n by : 
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We then calculate the ratio Ep^ of the current recording to the spectral model by the 



This calculation requires knowledge of M^^ whose calculation described below is 



10 carried out during cycle p-1. It comprises smoothing on the previous jfrequencies so as 
to limit the effect of the microseismic energy in its estimation. The proposed formula 
only uses the previous records and thus allows real-time conditions. In difficult cases, 
this sum can be replaced by an average or a current median. Deferred-time conditions 
are then necessary (the time required to carry out the second half of the recordings used 

15 in calculation of the average). 

If Mp+i^ represents the initial model of the active part of the recording n of cycle 
p+l, then : 



This model essentially consists of a weighted average of the current contribution 
20 and of the previous contributions. The formula used allows variation of the model : the 



relation : 
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contribution of cycle p-q to the model used for cycle p is multiplied by coefficient (1-h)** 
which tends to zero all the faster as h is close to 1. This allows model Mp^ to follow the 
seasonal variations, and ratios Ep^ to reflect only the daily variations. Knowing spectral 
model Mp.n and ratio Ep^ of the contribution to the previously established model, we 
5 deduce part Ap,n of the recording n of cycle p relevant to the active seismic part : 

Ap^ = Mp^ Ep^. 

Similarly, part Pp^ of the recording n of cycle p strictly relevant to the (passive) 
microseismic activity is calculated by the relation : 

10 The respective spectral contributions Cp.u of each source Si at the fimdamental 

frequencies and their harmonics, after completion of a cycle of measurements, provide, 
by inverse FFT, the seismograms (in the time domain) used within the context of active 
exploitation. 

The estimation of the underlying microseismic activity contained in record Rp^n is 
15 obtained by inverse FFT from the part Pp^ relevant to passive contributions. 

Updating of ratio Ep^ between the current contribution and the spectral model Mp^ is 
performed in a relatively fast way so as to allow fest variations of the signal. Updating 
of the initial model is performed in a much slower way; in fact, the initial model might 
not be updated. 

20 The interest of these updatings is to optimize estimation of the active signal when 

its variations cannot be controlled (seasonal and daily variations). In the most 
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favourable cases, they are not necessary and estimation of the active part amounts to the 
initial model. 

Another possible simplification consists in disregarding the seismic noise at the 
frequencies emitted (and, if need be, at the harmonics of these frequencies). One no 
longer undertakes to estimate the active part and one merely cancels or interpolates the 
lines corresponding to the frequencies emitted (and, if need be, those of their 
harmonics) in order to obtain the passive part. 

The harmonic lines of the frequency spectrum of each signal recorded are first 
suppressed, and the amplitude and the phase are interpolated at the line values at the 
frequency of the source; the record that has not been assigned by the emissions of the 
soiirces possibly containing a microseismic signal is then reconstructed by inversion. 
This implementation mode is rather suitable in the case of a limited number of sources 
and/or sources having a high signal variability witii time. 
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CLAIMS 

1) A method intended for active seismic monitoring of an underground formation 
(1) allowing to separate induced microseismicity signals from seismic signals emitted 
within the context of active seismic monitoring of an underground zone under 

5 development, comprising carrying out seismic recording cycles with emission of 
seismic waves in the formation by coupling therewith one or more seismic sources (5) 
which, in this case, emit simultaneously orthogonal signals so as to form a conq)Osite 
vibrational signal, reception of the signals reflected by the formation in response to the 
emission of seismic waves, recording of the signals received by at least one seismic 

10 pickup (2) and processing of the signals recorded so as to separate the respective 
contributions of the seismic sources to the signals received and to reconstract the 
seismograms equivalent to those that would be obtained by actuating the seismic 
sources separately, characterized in that the induced microseismicity signals are 
separated in the records from the seismic signals resulting from active monitoring 

15 operations, by isolating the contribution thereof by con5)arison with a reference spectral 
model by taking accovint of the spectral contributions of each source (Si) at the 
fundamental frequencies emitted and at the respective harmonics thereof and by 
reconstructing by inversion in the time domain the microseismicity signals. 

2) A method as claimed in claim 1, characterized in that the spectral contribution of 
20 the microseismicity signals to the spectrum of the signals received is obtained by 

subtracting the amplitude and phase values associated with the reference spectral model 
from the amplitude and the phase of the spectrum associated with the records. 
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3) A method as claimed in claim 1 or 2, characterized in that the reference spectral 
model is a current model formed by updating a previous spectral model by taking 
account of the spectral contribution of previous recording cycles. 

4) A method as claimed in any one of the previous claims, characterized in that the 
5 current spectral model is formed by determining a mean value of the frequency spectra 

formed from earlier and/or later records obtained for the same source and the same 
frequencies. 

5) A method as claimed in any one of the previous claims, characterized in that the 
current spectral model is formed by determining a median value of the frequency 

10 spectra formed from earUer records obtained for the same source and the same 
frequencies. 

6) A method as claimed in any one of the previous claims, characterized in that the 
current spectral model is formed by extrapolation or interpolation from the frequency 
spectrum from close spectral values. 

15 7) A method intended for active seismic monitoring of an underground formation 

comprising discrimination of the induced microseismicity signals from among signals 
emitted within the context of active seismic monitoring of an imdergroimd zone under 
development, comprising carrying out seismic recording cycles with emission of 
seismic waves in the formation by coupling therewith N seismic sources (Si) emitting 

20 simultaneously and controlled by orthogonal signals so as to form a composite 
vibrational signal, reception of the signals reflected by the formation in response to the 
emission of seismic waves, recording of the signals received by seismic reception 
. means'-and processing of the signals recorded so as to separate the respective 
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contributions of the seismic sources to the signals received and to reconstruct the 
seismograms equivalent to those that would be obtained by actuating the seismic 
sources separate^, characterized in that it conq>rises the following stages : 

a) for each recording n of a recording cycle p, tiie respective contributions (Cp,ui) of 
the various sources at the fundamental frequencies are calculated, 

b) the ratio (Ep^) of the contribution to a current spectral model (Mp,n) formed by 
updating a previous spectral model (Mp^.i) from fi^equencies emitted during the previous 
recording (n-1) and from their harmonics is then calculated, 

c) the part (Ap^ of tiiie recording n of cycle p that can be associated with the active 
seismic monitoring operations is deduced, 

d) the part (F^) of the recording n of cycle p strictly relevant to the passive 
microseismic activity is deduced, 

e) the seismograms that can be associated with die active seismic monitoring 
operations are formed by inversion in the time domain of the respective spectral 
contributions (Cpj.n) of each seismic source (Si) at the fimdamental frequencies and at 
their harmonics, after completion of a measuring cycle, and 

f) the underlying microseismic signals contained in the records are formed by 
inversion in the time domain from the part (Pp^) relevant to the passive microseismic 
activity. 

S) A method as claimed in claim 7, characterized in that the spectral contribution 
(Cp,i^) is obtained by multiplying a transfer function (Tp,vi,r) between a wavelet 
characteristic of the source and a seismogram associated with receiver r, by a wavelet 
characteristic of the source. 
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' 9) A method as claimed in claim 8, characterized in that said transfer fimction is 
continuously updated. 

10) A method as claimed in claim 9, characterized in that updating of said transfer 
function (Tp,,^) is obtained during a current cycle from an estimation (Tp,i^.i.,) made 
during a previous cycle and from an initial estimation (TOp^i^,) made during the current 
cycle by the relation : 

Tp,i^ = (1 - h)Tp^.,,r + hTOp.^n^. 



